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Abstract 



Electron motion in an oblique shock wave is studied by means of a one- 
dimensional, relativistic, electromagnetic, particle simulation code with full 
ion and electron dynamics. It is found that an oblique shock can produce 
electrons with ultra-relativistic energies; Lorentz factors with 7 ~ 100 have 
been observed in our simulations. The physical mechanisms for the reflection 
and acceleration are discussed, and the maximum energy is estimated. If 
the electron reflection occurs near the end of a large-amplitude pulse, those 
particles will then be trapped in the pulse and gain a great deal of energy. 
The theory predicts that the electron energies can become especially high at 
certain propagation angles. This is verified by the simulations. 
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I. INTRODUCTION 



Electron acceleration has been an important issue in plasma physics and astrophysics. In 
solar physics, for instance, the acceleration of electrons as well as of ions has received a great 
deal of attention; in solar flares, gamma rays with energies of several tens of MeV, which are 
emitted by the bremsstrahlung from the high-energy electrons, are often observed |l| , . In 
astrophysics, much more energetic electrons have been discussed. From the observations of 
X-rays and gamma rays, it is now believed that electrons with energies up to ~100 TeV are 
produced by a shell-type supernova remnant SN 1006 f^. It is also believed that in the 
Crab Nebula high-energy electrons accelerated up to ~100 TeV are produced by the pulsar 
wind . Theories based on the Fermi acceleration model have been proposed to account for 
these high-energy phenomena |^, 0]. As for laboratory plasmas, in an attempt to realize 
plasma-based accelerators, the electron acceleration in a relativistic space charge wave has 
been intensively studied by simulations and experiments M. 



Ion acceleration in a magnetosonic wave has been studied by many authors [^]- [|22 
In a single-ion-species plasma, the large positive potential formed in the shock region can 
reflect some ions and give great energies to them 1^- |18|. In a multi- ion-species plasma 
such as space plasmas, some of the majority ions (hydrogen) can also gain energies by 
the same mechanism. Furthermore, the transverse electric field in the shock wave can 



accelerate all particles of all kinds of (minority) heavy ions to nearly the same speed [jlQ 



In addition, in a turbulent plasma where particles can interact with many different 
large-amplitude magnetosonic pulses, some energetic ions can be further accelerated by 
a different mechanism These processes have been extensively studied by theory and 
particle simulation by many authors and have been applied to the production of high-energy 
particles in solar flares and interplanetary shocks 0- ||20[| . 

Strong electron acceleration in shock waves, however, has not been shown by particle 
simulations. In this paper we will theoretically and numerically study electron motion in 
a magnetosonic shock wave propagating obliquely to a magnetic field. We will show that 
an oblique shock can produce ultra-relativistic electrons. The preliminary result has been 
reported in Ref. [^. 

In Sec. 0, we outline some basic properties of oblique shocks. In Sec. |T|, we theoretically 
discuss the possibility of electron reflection. If electrons are reflected near the end of the 
main pulse, then after the reflection they can have great energies in the main pulse. 

Several authors discussed stationary pulse solutions for the nonlinear magnetosonic waves 
2^ - ||27|| . However, if the magnetic fleld is rather strong so that ujce/^pe ~ 1? where Uce 



is the electron cyclotron frequency and is the plasma frequency, or if the amplitude is 
quite large, then the pulse propagation is not perfectly stationary [|T^], [Q; small-amplitude 
fluctuations are generated in and around the main pulse. Even though the potential is 
positive in the main pulse in a shock wave, the electron reflection studied in this paper can 
occur in such a non-stationary pulse. 

In Sec. |V], we estimate the maximum energy of a reflected electron and discuss its 
dependence on plasma parameters. In Sec. |V|, we further study the electron acceleration by 
using a one-dimensional (one space coordinate and three velocity components), relativistic, 
electromagnetic, particle simulation code with full ion and electron dynamics. It will be 
shown that some electrons can be reflected near the end of the main pulse and are trapped 
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in the shock region. In the simulations, an appreciable number of electrons have energies 
7 > 100, where 7 is the Lorentz factor. Further, we will examine the dependence of 7 on the 
propagation angle and shock speed. It is verified that, as the theory predicts, there are special 
angles and shock speeds at which the electron acceleration is especially strong. Our work is 
summarized in Sec. |V^. In Appendix A, we describe relations among physical quantities in a 
magnetosonic wave. In Appendix B, we give a rough estimate for the magnitude of parallel 
electric field Eh on the basis of a simple physical picture. 



II. VELOCITIES AND FIELDS IN AN OBLIQUE WAVE 

We consider a magnetosonic wave propagating in the x direction with a speed Vsh in 
an external magnetic field in the (x, z) plane. We assume that the field quantities depend 
on X only; d/dy = d/dz = 0. Then, from the equation V ■ -B = 0, it follows that the 
X component of the magnetic field is constant, = B^q. The other components By and 
Bz are functions of x, and in the pulse region By as well as B^ can have finite values. In 
the wave frame where the time derivatives are zero [d/dt = 0), the electric field in the y 
direction is constant, Ey = EyQ, and Ez is zero. Hence we have B = {B^o, By, Bz) and 
E = {Ex, EyQ, 0). In Sections theoretical analyses will be made mainly in the wave 



frame. 

The velocity of the guiding-center position of an electron, Vg, may be written as 

Vg = {B/B)v\\ + Vd, (1) 

with Vd being the drift velocity 

Vd = cE X B/B^ -{cfi/e){B xVB)/B^. (2) 

Here c is the speed of light, — e is the electron charge (e > 0), v\\ is the velocity parallel to 
the magnetic field, and /x is the magnetic moment 

/i = mevl/i2B), (3) 

with me the electron mass and v± the gyration speed perpendicular to the magnetic field. 
In the far upstream region, we have no V-B-drift. Also, the z component of the velocity 
averaged over all the electrons in a small volume element must be zero, (vgzo) = 0; the 
subscript refers to the quantities in the far upstream region. Accordingly, from the z 
component of Eq. (|l|) we have the average parallel velocity 

CEyQ Bxo , . 

(^110) = -itlT^ 

in the wave frame (see Fig. |l]). Also, because {vgxo) = —Vsh, the x component of Eq. (|l]) 
gives 

Bxo I \ cEyoBzo ,_s 
Vsh = —^{v\\o} (5) 

-Do -Dq 

Thus EyQ is related to the shock speed Vgh through 
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Eyo = —VshBzo/ c. 



(6) 



For the definiteness, we assume that B^o, B^o and Vgh are all positive; and thus EyQ < 0. 

The nonlinear wave theory tells us the wave structure. A shock wave will have a positive 
electric potential (f{x). The quantities (p, B^, rii, and Ug have similar profiles 



see 



Appendix A); here n^, and Ue are the ion and electron densities, respectively. On the other 
hand, and By are proportional to the x derivatives of these quantities; for instance. 
Ex = —dip/dx. These relations among the quantities are obtained for small-amplitude 
waves. In the following theoretical analysis, we assume that the relations are also valid for 
large-amplitude waves. As we will see later in Fig. 0, simulation results also support this 
assumption. 

We denote by Xm the x position at which the potential takes its maximum value; hence, 
Ex{xm) — and By{xm) — 0. For later use, we also note that the quantity B^/B takes its 
maximum value at x = Xm', that is, if Bz{x) is smaller than Bz{xm), then 



Bla + BJxm? ^ Bl, + BJxY + BJxV ' 



^xQ^Bz{XmY B^Q + By 

This can be proved by the following equation 

Bz{xf _ [B,{xjj^ - Bz{xf]Bl^ + Bz{XmfBy{x) 



Bl^ + Bz{x^Y Bio + By{xf + Bz{xy [Bio + Bz{x^Y][Blo + By{xf + Bz^xf] 



(8) 



Because Bz{xm) is the maximum value of Bz{x), the term [Bz{xm)'^ ~ Bzix^jB^Q is positive. 
Hence the right-hand side of the above equation is positive. 



We know the stationary, finite-amplitude solutions for perpendicular waves |2J]- The 
stationary solitary wave solution is valid when the magnetic field is weak, Uce/^^pe << 1, and 
when the wave amplitude is not large, 1 < M < 2, where M is the Alfven Mach number; 
M = Vsh/vA with va the Alfven speed, v\ = Bq / {innimi) . Simulations and experiments 
show that large-amplitude pulses can propagate nearly steadily, even when the amplitudes 
are large or the magnetic fields are strong. Oblique pulses can also propagate nearly steadily. 
As the magnetic field becomes stronger or the wave amplitude is increased, non-stationary 
effects will become important |jl3|, 114 . 



III. ELECTRON REFLECTION 

Because the potential is usually positive in the shock region, some of the ions can be 
reflected when they go up the potential. As we will see later, electrons can also be reflected; 
its mechanism is, however, different from that of ion reflection. We give a physical picture 
of the electron reflection. 

From the nonrelativistic equation of motion for an electron 

we have an equation for the kinetic energy in the wave frame 
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nieiv'^ - vl)/2 = e{ip - ifo) - eEy^ j Vydt. (10) 

Here we have used the relation Ez = 0. The velocity v includes gyro, parallel, and drift 
motions. 

We define length s as 

ds = {B/B,o)vgxdt, (11) 
where Vgx is the x component of the electron velocity 

Vgx = {Bxo/B)v\\+Vdx, (12) 
(see Fig. ^. Thus ds/dt is given by 

ds B 

d^ = + b::'- ^''^ 

We can interpret ds as an infinitesimal length along the field line corresponding to the length 
dx = Vgxdt] ds = {B/BxQ)dx . In general the guiding-center velocity Vg is not parallel to 
the magnetic field B. Hence the guiding center does not move along the field line. For a 
stationary, one-dimensional problem {d/dt = and d/dy = d/dz = 0), all quantities depend 
only on x. As a result, for any stationary function f{x) we have 



/ f{x){B/Bxo)dx= / /(s)ds, (14) 



where Si and S2 are the positions along the same field line; their x positions are Xi and X2, 
respectively. 

If we define quantity F as 

E|| = -dF/ds, (15) 
then with the aid of (|llj) and the relation E\\ = E ■ B/B, we have 

F = ~ f (M^!^±3M^dx, (16) 

which can be written as 

Using the electric potential and the vector potential A [B = V x A), we find 

F = ^-{vsh/c){Bzo/Bxo)Az, (18) 
where Az is the z component of the vector potential, 

Az = - [ Bydx. (19) 
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The quantities if and have similar profiles, while By is proportional to dB^jdx p8[, p9 . 
Hence and F have profiles similar to ^p. 

If we define quantity F; in the laboratory frame (the subscript / denotes the laboratory 
frame) as 

F; = -| Ei\\{BilBi.,^)^xu (20) 

then for a stationary, one-dimensional problem, that is, for a case where wave profiles can 
be written as f{xi,ti) = f{xi — Vshti), we find 

F. = -/(E..-^^^^§^)dx.. (21) 

One can readily show that it is related to F in the wave frame as 

F = ishFi, (22) 

where 'jsh is defined as 

lsH = {l-v'JcY'^'. (23) 



In Appendix B, we give a rough estimate of the maximum value of F, Eq. (|B9|) , found from 
qualitative physical considerations. 

Combining Eqs. ( p!OD and (p!8|), we eliminate (f 

met;V2 - = e{F - Fq) - {eEyo/ B^o){A, - A,o + J VyB^odt). (24) 

By virtue of the definition of A^, (0), and the relation dx = Vxdt, this equation can be 
written as 

mevy2 - myj2 = e{F - Fq) + e{Eyo/B.,o) J {v,By - VyB,o)dt. (25) 
Then, using the z component of the equation of motion, we find 

- mevl/2 = e{F - Fq) - mec{EyQ/ Bx,q){vz - v^o)- (26) 
In the drift approximation, the kinetic energy can be expressed as 

met;V2 = me(t;J + i;^)/2 + /i5. (27) 
Hence, averaging Eq. (^) over the electron cyclotron period, we have 

{me/2){vf\ + vl) + fiB + mec{Eyo/Bxo)vgz - eF 
= (me/2) (f Jo + vIq) + fxBo + mec{Eyo/Bxo)v^o - cFq. (28) 

With the aid of the z component of Eq. (|T]), we can eliminate Vgz and obtain the following 
equation 
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{mj2){v\\ - Vrvf + K = e{F- Fq) - fi{B - Bq) 
+ (me/2)(t;Jo + vIq) - m^c{Eyo/ B.^o){vdz - v^o) - inie/2)vj, (29) 

where Vrv and K are defined as 

= -cEyoBj{B,oB) (> 0), (30) 



K = -my^j2 (<0). (31) 
If we eliminate EyQ by substituting Eq. (^ in Eq. (|30D , be expressed as 

Vrv = VshBzoBz/{B^oB). (32) 

We will see below tliat tlie velocity Vgx is reversed wlien v\\ = Vrv In the pulse region, the 
quantity e{F — Fq) can be much larger in magnitude than the other terms on the right-hand 
side of Eq. (|29| ) (see Appendix B). We denote the left-hand side of Eq. (^) by Eg, i.e., 
Eg = (me/2) (f II — Vrv)'^ + K. As shown in Fig. |], at a fixed point x, Ee has its minimum 
value K when t;|| = Vrv It is negative in the region 

< < 2Vrv (33) 

Because Vdx = cEyoBz/B^ , it follows from Eq. ([I2|) that Vgx = if v\\ = Vrv As can be 
seen from Eq. (^), most of the electrons have negative parallel velocities in the far upstream 
region. Accordingly, we suppose that the initial parallel velocity is negative, v\\ < 0. If 
f II gradually increases and exceeds Vrv at a certain position, then the x component of the 
guiding center velocity, Vg^, changes from negative to positive values there, i.e., it is reversed. 

As an electron moves, i.e., as x changes, the quantities K, F, B, etc. also vary. At any 
point X, K is the minimum value of the left-hand side of Eq. (^). If there is a region where 
the values of the right-hand side of Eq. ( PP| ) become smaller than K, then the electron 
cannot enter there; it will be reflected. We recall that ip, F, and Bz have similar profiles 
and that at x = Xm the quantity (Bz/B) has a maximum value [see Eq. (^]. Because K is 
proportional to {Bz/BY, it will be quite small {\K\ is large) around the point x = Xm, and 
F is large. In front of or behind the main pulse, F is relatively small and K is relatively 
large; e{F — Fq) could be smaller than K. Thus, if the electron reflection occurs, it will be 
in such regions. 

In quasi-perpendicular shocks, where (Bzo/B^o) » 1, Vrv is large. Hence K has large 
negative values. Thus the electron reflection in this mechanism will be difficult. 

Figure ^ shows a schematic diagram of the trajectory of guiding center. Here it is assumed 
that the reflection takes place at point D; the dotted line shows the orbit of a passing electron. 
As an electron moves from point A to C, it moves to the negative y direction because of the 
Ex X Bz drift. It gains kinetic energy AEi from the electric potential 

AE, = eifiixc) - eifixA) (> 0). (34) 

At the same time, it loses energy AE2 because of the electric field Eyo 

AE2 = -eEyo{yc-yA) (< 0). (35) 
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The net change in the energy is therefore 



AE = AEi + AE2. (36) 

Even though the magnitudes of AEi and AE2 are quite large, they almost cancel when an 
electron moves with drift approximation; in particular, in a perpendicular pulse they have 
exactly the same magnitude and AE = |^ . However, if an electron is reflected and moves 



from D to E, then it would gain energies from both E^ and EyQ. As a result, the increase in 
energy is 

AE = AEi + AE3, (37) 

where Ai^s is defined as 

= -eEyoivE - Va) (> 0). (38) 



IV. ESTIMATE OF ENERGY INCREASE 

Next, we will obtain the maximum energy of a reflected electron. Because it can have 
quite high energy, we will use the relativistic equation of motion 

where 7 is the Lorentz factor. For the stationary, one-dimensional system being considered, 
this can be integrated to give the energy conservation equation 



me7C^ - me7oC^ 



e(v? - v^o) - eEyo J Vydt. (40) 



This is a relativistic form of Eq. (IC). Again, using the quantities F, Az, and the equation 



of motion, we can rewrite Eq. ( pH ) as 

me7c^ - me7oc^ = e(F - Fq) - mec{Eyo/ Bm){iVz - JqVzo). (41) 
If we define quantity h as 

h = nieC^ - nieVshBzovJ Bm, (42) 

then from Eq. (|4l|) we have 

7= [e(F-Fo) + /io7o]//i. (43) 

Equation (|^) indicates that h is positive if Bzo/B^o is of order unity. If Bzo/B^o is much 
greater than unity, then h can be negative. In addition, Eq. (|43|) shows that, if the initial 
value of h is positive, > 0, then h is always positive in the region where {F — Fq) is 
positive. 
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We show in Fig. ^ schematic diagram of the function [e(F — Fq) + /io7o]- We suppose 
that F has a maximum value at x = Xm and has a minimum value right behind the main 
pulse. If the reflection occurs, it would be in this dip, as suggested in the previous section. 
The reflection point will be denoted by Xr- In the top panel this function e(F — Fq) + /io7o is 
always positive. In the second one it becomes negative right behind the main pulse; however 
its maximum value is positive {ho can be either positive or negative). In the bottom one 
it is always negative. For the present acceleration mechanism, the top panel is the most 
important; the strong acceleration takes place in this case. 

Such dips would not be present in the perfectly stationary solutions EM- E^. In large- 



amplitude magnetosonic waves, however, non-stationary dips can be generated |T^, [|14 
We here just assume that there is a dip behind the main pulse and discuss its effects on 
particle orbits. 

First we consider the top panel in Fig. |^: An electron with positive ho is reflected at a 
certain point near the end of the main pulse. We will show that, after the reflection, h has 
its minimum (positive) value at x = Xm', hence 7 has its maximum value there. 

As mentioned earlier, the parallel velocity v\\ is already positive when an electron is 
reflected, i.e., when Vgx is reversed, even if its initial value is negative. After the reflection, 
f II will be increased by the (positive) parallel electric force. Because E\\ changes its sign at 
X = Xm, v\\ will have a maximum value ai x = Xm (unless changes to negative values again 
before reaching the point x^ ) • The parallel velocity will become the dominant component in 
the electron velocity, v v\\B/ B. The quantity Bz/ B as well as ip and Bz has a maximum 
value at this point. Thus Vz-, which can be approximated as Vz ~ v\\Bz/B, will also have the 
maximum value there. Because ho and [Fm — Fq) are both positive, h is positive at x = Xm- 
In addition, Eq. (|4^) indicates that h decreases (draws closer to zero) as Vz increases. Hence, 
h of a reflected electron will have a minimum positive value at x = x^. Its energy 7 will 
therefore have a maximum value there. If h can be close to zero, then 7 would become quite 
large. 

If Bxo is small and Vzo has a large positive value, then ho becomes small (can be even 
negative). We now consider such particles; the second panel in Fig. ^ Also in this case, 
Vz will have larger values after the reflection; thus h will decrease. At the reflection point 
in the dip, the quantity [e(F — Fq) + /io7o] is negative. Consequently, h{xr) < 0. If the 
reflection occurs, the electron would move in the positive x direction. Because Vz increases, 
h will further decrease. Hence, the particle would never be able to penetrate the region 
again where /17 > 0. The strong acceleration in this mechanism, therefore, will not take 
place. 

Next, we consider a case where the quantity [e{F — Fo) + /io7o] is always negative (the 
bottom panel in Fig. Then, h is always negative. If the reflection occurs at x = x^, 
then \h\ would increase while |e(F — Fq) -|- /io7o| would decrease. Hence, even if the particle 
reaches the point x^, 7 would not have an extremely large value. Therefore, the strong 
acceleration is not expected. 

We can estimate the maximum value 7^ in the wave frame by substituting the maximum 
value of F, Eq. (|B12| ), in Eq. (^31). The maximum 7 in the laboratory frame can be obtained 
by the Lorentz transformation as 

llm = 1sh{l + VccVsh/ C^hm - Ishim- (44) 
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Now let us consider a case where h becomes quite small, /i ~ 0, at x = Xm- There the 
electron speed would be close to the speed of light c. Hence, roughly, we have ~ cB^/B. 
The quantity h can then be expressed as 

h = nieC^ ^1 ~ tan 6*^ , (45) 

where tan 6* = Bizq/Bi^q. Here we have used the relations B^o = 'jshBizo and B^q = Bi^o- 
Because B^ becomes quite large in the shock region, we may put as Bz/B ~ 1. Equation 
(^) then suggests that the electron energy can become extremely large when 

Jsh{vsh/c)tan9 = 1, (46) 

or, equivalent ly 

c cos 9 = Vsh. (47) 



V. SIMULATION STUDIES 

A. Simulation Method 

We further study the electron acceleration by using a one-dimensional (one space co- 
ordinate and three velocity components), relativistic, electromagnetic, particle simulation 
code with full ion and electron dynamics |2T|. The total system length is = 4096 



for most of the simulation runs, where Ag is the grid spacing. All lengths and velocities in 
the simulations were normalized to Ag and UpeAg, respectively, where Upe is the spatially 
averaged plasma frequency. 

We use a bounded plasma model. The plasma is limited to the region 400 < x < 3696. 
The particles are specularly reflected at x = 400 and at x = 3696. The radiation leaving the 
plasma region is absorbed in the vacuum regions, < x < 400 and 3696 < x < 4096. Thus 
the electromagnetic interactions between the two plasma boundaries through the vacuum 



regions are made negligibly small |T2 



In most of the plasma region, i.e., in 500 < x < 3696, each particle species initially has 
a uniform density and has a Maxwellian velocity distribution function. In the small region, 
400 < X < 480, the plasma density is four times as high as that in the main region. That 
is, the initial plasma density can be written as 

He = 4no for 400 < x < Xi, 

Ue = 4noexp[— (x — xi)'^/{2(f )] for xi < x < 500, 

He = no for 500 < x < 3696. 

Here xi = 480 and d = 12. In the high-density region, 400 < x < 480, particles have 
shifted Maxwellian velocity distribution functions with average velocity Vp^; fj ~ exp[— (v — 
Vps)^ / {2v^j)], where the subscript j refers to particle species. These particles act as a 
piston. They push the neighboring particles and excite a shock wave. We can change the 
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shock strength by changing the magnitude of Vps] the velocity Vps is perpendicular to the 
external magnetic field so that no particles have extremely large parallel speeds initially. 
More detailed description about the simulation code can be found in Ref. . 

The simulation parameters are as follows. The number of simulation particles is Ni = 
Ne = 262,144. The ion-to-electron mass ratio is mi/me = 100. The ratio ujce/^pe is 
iOcej^pe = 3.0 in the far upstream region. The electron and ion thermal velocities are 
VTe/c = 0.38 and vti/c = 0.01, respectively. The Alfven speed is va/c = 0.3. The electron 
skin depth is c/upe = 4Ag. The time step is sufficiently small, WpeAt = 0.02, so that At is 
much smaller than the plasma and cyclotron periods even in the shock region. As in the 
theoretical model, the external magnetic field is in the {x, z) plane, and waves propagate 
in the x direction. 



B. Simulation Results 



First, we study a shock with the propagation angle = 45°; i.e., B^q = B^o- Figure |^ 
shows profiles of Bz in the shock wave at various times. Here, the magnetic field is normalized 
to BzQ. The Mach number is observed to be M = 2.3, i.e., the shock speed is Vsh = 2.3fA. 
The maximum value of Bz is 4~5 times as large as Bzq. The main pulse propagates nearly 
steadily. However, the profile is not perfectly stationary. 

Figure |^ displays field profiles at ujpet = 680; here the magnetic and electric fields are 

], Ey and Bz have similar profiles, 
while Ex, Ez, and By are quite small in magnitude at the point where Bz has its maximum 
value. Figure ^ shows phase space plots of electrons at Upet = 680. In the shock region, 
ultra-relativistic electrons are produced. The maximum value of Lorentz factors is 7 ~ 130. 
We also see that and Pz have greater values than py, which is consistent with the picture 
that the parallel velocity is dominant in the motion of high-energy electrons. From Figs. |^ 
and D, we see that some electrons are refiected near the end of the main pulse and that the 
high-energy electrons are present in the shock region. 

Figure ^ shows profiles of -Ey, F, and the electric potential ip at Upet = 650 and at 
ujpet = 680. The quantities F and (p denote eF/{meC^) and eip/{meC^), respectively. The 
quantities F and were obtained in the simulation by the following equations 

F{x) = -fE\\{B/Bxo)dx, (48) 



both normalized to Bq. As the theory predicts [Bq], |29 



^(x) = - / Exdx. (49) 



As mentioned earlier, roughly speaking, the parallel electric field E\\ is positive in the region 
where dBz/dx and dF/dx are negative. It is negative in the region where they are positive. 
Comparing the plots at different times, we see that the wave profiles are not perfectly 
stationary. They vary with time with small-amplitudes. In particular, as shown in the plot 
at ujpet = 680, F sometimes becomes negative near the end of the main pulse, where the 
electron reflection can take place. 

Next, we will study trajectories of electrons; we will find that it is the refiected electrons 



that gain great energies. Figure |T0] displays time variations of {x — Vght), y, z, and 7 of 
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electrons in a shock with Vsh = I-TSiia; {x — Vsht) is the x position in the wave frame. The 
lengths are normalized to the electron skin depth c/upe. (Here the system size is a half of 
that in the previous run; the other parameters such as N^/L^ are the same.) Two typical 
examples are presented here. The thick lines show an electron that was accelerated, and 
the thin lines represent an electron that was not accelerated; even the low-energy electron 
here has energy 7 ~ 10 in the shock region and 7^5 behind the shock. In all of the four 
panels, the low- and high-energy electrons both have small-amplitude, cyclotron oscillations 
with short periods, Upf^t ~ 20; note that the cyclotron period depends on 7 as well as the 
field strength B. At the same time, the accelerated electron has a large-amplitude, long- 
period oscillation with period Up^t ~ 170; it begins when the particle enters the shock, 
ujpet = 180, and continues until the end of this simulation run. (Simulations show that 
this period increases with 7. To observe this oscillation as many times as possible, we have 
chosen here the case where 7 is not so large.) For comparison with Fig. ^ we show the 
times corresponding to points B, C, D, and E. From the long-period oscillation of (x — Vsht) 
in the top panel, we see that the electron that will become high energy is reflected and then 
trapped by the shock wave. On the other hand, the quantity {x — Vsht) of the low-energy 
electron keeps decreasing, which means that this particle passes through the shock region 
without strong interactions. The second panel shows that the y position of the accelerated 
electron also oscillates, while the third one shows that z increases rapidly after the reflection. 
The average value of velocity over the period from ujpet = 220 to 400 is v^/c ^ 0.6, and 
the value of Vz at maximum 7 is Vz/c ~ 0.8. The maximum 7 in the bottom panel is 7 ~ 40. 

By the Lorentz transformation of these data, we can find v and B in the wave frame. 
At the time and position of maximum 7, they are v = c(0.27, 0.13, 0.95) and B = 
i?o(0.70, 0.16, 2.56). This confirms our statement in the previous section that is the 
dominant component in v there; i.e., Vz — v\\Bz/ B. 

According to the theory, (0), electron energies can be extremely high at some angles 
and shock speeds. To examine this prediction, we have carried out several simulations with 
different values of the propagation angle 9, keeping other parameters unchanged; Lr^ = 
4096Ac,, Uce/ujpe = 3.0, and Vsh — 2.1va- Figure |Tl| shows the maximum value of 7 as a 
function of the propagation angle 6. Here the electron energy takes the highest value at the 
angle 6 ~ 52°. For these simulation parameters, Eq. (^Tf ) predicts that 7 takes a peak value 
at ^ ~ 53°. Thus the theory and simulation are in good agreement. Further, we carried out 
simulations with different values of the shock speed Vsh', the propagation angle is fixed to be 
6 = 45°. Figure O shows the maximum value of 7 as a function of the shock speed Vsh- The 
observed values have a peak at Vsh/c — 0.7. For these parameters, Eq. (|^) predicts that 7 
takes large values for Vsh/c ^ 0.71. Again, the theory and simulation are in good agreement. 



VI. SUMMARY 



We have studied electron motion in a shock wave propagating obliquely to a magnetic 
field. First, we analytically discussed electron motion in an oblique shock wave. It is 
pointed out that, if an electron is reflected near the end of a large-amplitude pulse, it would 
gain a great deal of energy from the electric field formed in the wave. The condition for 
the reflection is then examined. Further, the maximum energy of a reflected electron is 
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estimated, and its dependence on plasma parameters is discussed. Next, we investigated the 
shock evolution and associated electron acceleration by using a one-dimensional, relativistic, 
electromagnetic, particle simulation code with full ion and electron dynamics. It was shown 
that an oblique shock can produce ultra-relativistic electrons; Lorentz factors with 7 ~ 100 
have been observed. As the theory predicts, at certain propagation angles electron energies 
become extremely high. The electron reflection and resultant acceleration takes place when 
the wave amplitude is large and the magnetic field is rather strong, Uce ~ ^pe- 

In the present paper, strong electron acceleration has been demonstrated. As future 
work, it would be desirable to develop further the quantitative theory for large-amplitude 
oblique waves. The evaluation of the electric field strength along the magnetic field will be 
especially important, which will enable us to estimate the maximum energy more accurately. 
From the view point of nonlinear wave propagation, motion of trapped electrons and their 
effects on the wave evolution are also quite interesting. 



APPENDIX A: RELATIONS AMONG QUANTITIES 

We here describe relations among the physical variables, using the approximation appro- 
priate for low-frequency magnetosonic waves with small amplitudes (for the details of the 
calculations, see Ref. |2^). We assume quasi neutrality, n ~ ~ ne and use the following 
stretched variables 

r = e^/^t, (Al) 

i = e^'\x-v,ot). (A2) 

(fpo is the wave propagation speed in the long-wavelength limit in a finite beta plasma) and 
expansion 

n = no + eni + e^n2 + (A3) 

Ey = eEy^ + e^Ey2 + (A4) 

B, = B,o + eB,i + e^B,2 + ..., (A5) 
then, we can express the lowest order perturbations in terms of ni as 

Bzi/Bo = cEyi/{BoVpo) = [{vL - c1)/vlsm9]{ni/no), (A6) 



cE^i _ VpoByi _ {uce - tOci) (^po - cj,)vlo cot e dni 
Bo Bo uJceUJci (fpo - cos^ 6) no 



(A7) 



Here Cg is the sound speed, = (FePeo + ^iPio)/ {nonii) with Tj (j = e or z) the specific heat 
ratio and pjo the equilibrium pressure. We can obtain KdV equation by proceeding to the 
next order, 0(e^/^). 

The relations ( [A6| ) and ( |A7| ) are valid, even when the wave profile is not soliton-like. 
That is, perturbations propagating in the same direction with Vx ~ fpo have the relations 
and (1X71). 
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APPENDIX B: MAGNITUDES OF En AND F 



For a nonrelativistic, small-amplitude oblique magnetosonic wave, the electric field par- 



allel to the magnetic field, E\\, has been analytically obtained [Q. However, for a large- 
amplitude wave, it will be quite difficult to obtain E\\ in a rigorous manner. Here, on the 
basis of a simple physical picture, we will give a rough estimate of the parallel electric field 
E\\ in a relativistic, large-amplitude, oblique shock wave. Further, by using it, we will obtain 
the maximum value of the quantity F. 

We consider in the wave frame. In the fiuid model, the equation of motion may be 
written as 



m 



'J 



dt 



qjE + q,^ . (Bl) 



Here the pressure term is neglected. The subscript j refers to ions j = i ot electrons j = e. 
The time derivative d/dt is defined as 



We denote by b the unit vector along the magnetic field. Then multiplying Eq. ( |B1|) by b, 
we have 



m 



d(7j6 • Vj) db 



■J 



^ni'"r-n = (liE\\. (B3) 



dt ' dt 



Let ti be the time when a fiuid element is at the leading edge of the shock (the location 
where the shock profile begins to sharply rise) and ^2 be the time when it is at the point 
Xm'-i Xm is the location where the electric potential and magnetic field take their maximum 
values. Then, integrating from ti to t2, we find for the ions 

A 

rni(7i2&2 ■ - 7ii&i • ^'ii) - mi{jiVi) ■ (62 - bi) ~ eE\\- — -. (B4) 

\{Vix}\ 

Here the bracket indicates the mean value over the time period (^2 — ^1), and A is the 
shock width. We have a relation (t2 — ^i) — ^/\{vix)\- For oblique shocks, A will be of 
the order of the ion inertial length, A ^ c/ujpi. The velocity at time ti is vn = —Vgh^x, 
where is the unit vector in the x direction. The velocity Vix will be slowed down by 
the longitudinal electric field. Thus Vix2 will be smaller in magnitude than Vixi- Because 
Bz is large at the point Xm and By{xm) ~ 0, the unit vector 62 is nearly parallel to the z 
direction. Consequently, the inner product (62 ■ ^^{2) will be quite small; here we neglect the 
z component of Vi2. Since the x component of 61 is given by B^o/Bq, we have 

friii'jnVsh + {7iVix))Bxo/Bo ~ eEy A/| (wj^) |. (B5) 

This gives the strength of the parallel electric field as 



A Vsh \ Vsh J Bq 
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Because \vix\ is smaller than Vsh, the range of \{vix)\ is expected to be 1/2 ~ \{vix)\/vsh < 
1, except maybe for extremely strong shock waves. In fact, if we assume a simple time 
dependence of Vix 

Vix = -Vsh[l - ai{t - ti)/(t2 - ti)] for ti<t< t2, (B7) 

with ai a constant (0 < Oi < 1), then, after integrating over time, we have {vix) = —Vghi^ — 
ai/2). 

The Lorentz factor for the fluid ion will be close to unity. Accordingly, for instance, if 
the ion fluid velocity Vix is decreased to a value ~ —Vsh/'2, then the value of the content in 
the parentheses would be about 1/4. 

The electron speed will be ~ Vgh when the electrons encounter the shock and will be close 
to the speed of light when they are at the point Xm- If we assume that the electron velocity 
is nearly parallel to the field line in the shock region, then [v ■ db/dt) can be neglected. 
Hence we have 



me{-'ye2C + JelVshBxo/ Bq) ~ -eE|| (ts -ti). (B8) 

The second term on the left-hand side, 'jeiVshBxo/ Bq, can be neglected compared with the 
first term 7e2C. The ions and electrons will have nearly the same time period (t2 ~ ^i) to 
keep charge neutrality. Substituting Eq. (BB) in ([B8|) yields 



7e2 ^ 7ii + (B9) 

rrieC \ Vsh ) Bq 

At be estimated as Vz — —cBz/B (for reflected electrons the sign is reversed). 

For a large-amplitude shock, the ratio Bz/B is an order-unity quantity there. The quantity 
h defined by Eq. (|42|) is obtained as 



h2 ~ mec2[l + VshBzo/icBxo)]. (BIO) 

If we average Eq. (^) over electrons in a small volume element, we have 

{h^) = e{F-Fo)+meC^^o). (Bll) 

Here we used the relation {vzo) = 0. If electrons have the velocity Vz — —cBz/B at x = Xm, 
then we find the maximum value of F from Eq. (^) as 

eiFra - Fo) ^ /l27e2 " m,c^{^o). (B12) 

In many practical cases, we can take (70) to be order unity. 

Roughly speaking, the quantity e(Fm — Fq) is of the order of rriiCVsh- It is, therefore, 
much greater than electron thermal energy rUeV^^ (~ fiB) and drift kinetic energy rngV^ 
(~ m^v^^fj. 
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FIGURES 



FIG. 1. Velocities and fields in the far upstream region in the wave frame. The velocity Vd is 

Eyo X Bq drift, and —v^h is equal to average guiding-center velocity in the x direction {vgxo)- 
FIG. 2. Schematic diagram of magnetic field, electron velocity, dx, and length ds. 
FIG. 3. Plot of function E(,{v\\) = (me/2)(t;|| — Vrv)"^ + K at a fixed x position. 
FIG. 4. Schematic diagram of electron orbit in (x, y) plane. 
FIG. 5. Schematic diagram of quantity [e{F — Fq) + /iqTo]- In the top panel, this quantity 
is always positive. In the second one, it is negative in the dip. In the bottom one, it is always 
negative. 

FIG. 6. Magnetic field profiles of an oblique shock at various times. 
FIG. 7. Snapshots of field profiles. Electric and magnetic field profiles at Upet = 680 are plotted. 

They are normalized to Bq. 

FIG. 8. Phase space plots of electrons. 
FIG. 9. Profiles of £"11, F, and (p at two different times. Here, F and (p denote eF/{meC^) and 
e^p / {m,(.c^) , respectively. 

FIG. 10. Time variations of {x — Vght ), y, z, and 7 of electrons. The lengths are normalized to 

the electron skin depth c/upe- 

FIG. 11. Maximum 7 versus propagation angle 9. The propagation speeds for these shocks are 

Vsh ^ 2.lvA- 

FIG. 12. Maximum 7 versus shock speed Vgh- The propagation angle 6 is fixed to be ^ = 45° 
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